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ABSTRACT 



A torus around a stellar mass Kerr black hole can emit about 10% of 
the spin-energy of the black hole in gravitational radiation, potentially asso- 
ciated with a gamma-ray burst. Wide tori may develop buckling modes by 
the Papaloizou-Pringle instability and gravitational radiation-reaction forces in 
the Burke-Thorne approximation. Gravitational wave experiments may discover 
CN) . these emissions in a fraction of nearby supernovae. This provides a test for Kerr 

black holes, and for GRB inner engines by comparison with the de-redshifted 

CN . durations of long GRBs. 
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1. Introduction 



Stellar mass black holes surrounded by a compact torus may represent catastrophic 
events such as core-collapse in a massive stars and black hole-neutron star coalescence. 
These scenarios have been considered as sources of cosmological gamma-ray bursts (Woosley 
1993; Paczynski 1991). We may consider black hole-torus systems and their emissions more 
generally, especially when the black hole is rapidly spinning. Their potential association 
with GRBs provides observational constraints on their evolution. 

A torus around a rapidly rotating black hole may develop a state of suspended accretion 
for the lifetime of rapid spin of the black hole (van Putten & Ostriker 2001). This points 
towards major energetic output in "unseen" emissions, in gravitational radiation, magnetic 
winds, thermal emissions and neutrino emissions (van Putten & Levinson 2002). The energy 
E gw emitted in gravitational radiation is expected to be about 10%E rot , i.e., 

E gw ~ 6 x 10 53 erg (1) 
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for a 1OM black hole. This output (1) may be detected by the upcoming Laser Interferomet- 
ric Gravitational Wave Observatory LIGO (Abramovici et al. 1992) and the French-Italian 
counter part VIRGO (Bradaschia et al. 1992), possibly in combination with any of the bar 
or sphere detectors currently being developed. This provides a calorimetric compactness test 
for Kerr black holes (van Putten 2001), and a means of identifying the inner engine of GRBs 
by comparison with de-redshifted durations of long GRBs. 

In the Woosley-Paczynski-Brown scenario of hypernovae (Woosley 1993; Paczyhski 1998; 
Brown et al. 2000; Lee et al. 2002), core-collapse in rotating massive stars forms a Kerr black 
hole surrounded by a compact disk or torus. Long GRBs correlate with star-forming regions 
(Bloom et al. 2000) and, hence, young massive stars, possibly in binaries. The rotating 
black hole may produce wide-angle ejecta back into the interstellar medium leaving behind 
a soft X-ray transient with a chemically enhanced companion star (Brown et al. 2000), such 
as GRO J1655-40 (Israelian et al. 1999) and V4641Sgr (Orosz et al. 2001). This potential 
supernova and SXT association is important in identifying progenitors to GRBs and their 
inner engines. The beamed output of true GRB energies of E 1 = 10 50_51 ergs (Frail et al. 
2001) represents a minor energetic output for a long GRB from a rotating black hole. The 
potential for long gamma-ray bursts from rotating black holes suggests that GRB associated 
supernovae may emit bursts of gravitational radiation, e.g., GRB 980425/1998bw (Galama 
et al. 1998), GRB 011121 (Bloom et al. 2002) and GRB 011211 (Reeves et al. 2002). 

In this Letter, we suggest LIGO/VIRGO searches for bursts of gravitational radiation 
from back hole-torus systems using upcoming continuous all-sky supernovae surveys. Focus- 
ing on supernovae may serve to reduce data analysis by their well-determined coordinates, 
which includes distances. The expected gravitational wave-spectrum is here identified with 
multipole moments in a wide torus due to a Papaloizou-Pringle instability, by extension 
of the theory for slender tori (Papaloizou & Pringle 1984; Goldreich et al. 1986) and in- 
cluding the secular effect of gravitational radiation backreaction-forces in the Burke-Thorne 
approximation (Thorne 1969). 

2. GRB/GWBs from rotating black holes 

A rapidly rotating black hole surrounded by a magnetized torus may develop a state 
of suspended accretion (van Putten & Ostriker 2001). An equilibrium magnetic moment in 
its lowest energy introduces a high incidence of the black hole luminosity into the inner face 
of the torus through a surrounding magnetosphere with slip/no-slip boundary conditions, 
and enables support to an open flux tube to infinity with slip/slip boundary conditions (van 
Putten & Levinson 2002). 
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Fig. 1. — Shown is the histogram of redshift corrected durations of 27 long bursts with 
individually determined redshifts from their afterglow emissions (sample from Djorgovski et 
al., 2001, astro-ph/0107535, and references therein, and updated (Djorgovski, 2002; Hurley, 
2002). The mean value of the durations of all bursts is 38s; the average is 23s without the 
long bursts GRB 980703 (T 90 /l + z = 400s) and GRB 000911 (T 90 /l + z = 243s), which 
consist of two (in BATSE; one in Ulysses), respectively, three well-separated sub-bursts. 
GRBs from rotating black holes are expected to be accompanied by "unseen" emissions in 
gravitational radiation from the torus. This predicts a similar distribution of durations for 
their bursts of gravitational radiation with an expectation value of about one-half minute. 
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Differential frame-dragging on an open flux tube may create a minor energy output Ej 
in low-sigma baryon-poor jets. In the proposed association to GRBs, Ej = E 1 /e represents 
a few permille of the rotational energy of the black hole, where e ~ 0.15 denotes a fiducial 
efficiency of the conversion of kinetic energy-to-gamma rays. The associated horizon half- 
opening angle 9h — 35° provides a bound for the observed half-opening angle 9j on the 
celestial sphere, following collimation. A major energy output is expected to derive from 
the torus, upon catalyzing spin-energy from the black hole into gravitational radiation (1), 
winds, thermal and MeV neutrino emissions. We next consider the formation of a multipole 
mass moments in the torus, which defines its spectrum of gravitational radiation. 



3. Multipole moments in wide tori 

The effect of shear on the stability of a three-dimensional torus of incompressible fluid 
can be studied about an unperturbed angular velocity Q = Q a (a/r) q , where q > 3/2 denotes 
the rotation index. In the inviscid limit, we may consider, by Kelvin's theorem, irrotational 
modes in response to irrotational initial perturbations to the underlying flow (vortical if 
q 7^ 2). We expand their harmonic velocity potential 

= £„a„(r,M)z n , A0 = O, (2) 

in cylindrical coordinates (r, 9, z). The equations of motion can be conveniently expressed in 
a local cartesian frame (x, y, z) with angular velocity Q a = M 1 l 2 ar z l 2 at r = a about a central 
mass M, where x = r—a, d x = d r and d y = r~ l do. (We can readily switch between cylindrical 
and cartesian coordinates in coordinate invariant expressions.) Infinitesimal perturbations 
oc e ime ~ lU},t of frequency uo' as seen in the corotating frame at r = a satisfy the linearized 
equations of momentum balance. For an azimuthal quantum number m and on the equatorial 
plane z — 0, they are, in the notation of (Goldreich et al. 1986), 

-iau-2Qv = —d r (h + &), , , 

-iav + 2Bu =-ik(h + <5>), ^ 

where h denotes a perturbation of the unperturbed enthalpy d r h e = Vt 2 r — M/r 2 about 
z = 0, $ a potential, IB = (2 — q)il, k = m/r, and a = uo' — m5VL (5VL = VL — Q a ) produces 
the Lagrangian derivative D t = —ia + ud x . The equation of motion for the vertical velocity 
component, —iaw = —d z h, is decoupled in z = 0. In earlier linearized treatments (Goldreich 
et al. 1986), variations of 2B across the torus are neglected, limiting the discussion to narrow 
tori defined by h e (x±) = 0. For wide tori, we here include (2B) X = —(q/r)2B. The equations 
of motion (3) obtain 

{d 2 r + r- x d r - m 2 r~ 2 )a = qr~ l a' (4) 
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for an azimuthal mode number m. Solutions symmetric about the equatorial plane satisfy 

2 

<P = a -^-d r a + O(z 4 ), a = r p + + \r p - , (5) 
2r 

where p± = q/2± ^/q 2 /^ + m? and A is a constant. 

The unperturbed inner and outer boundaries x — x± of the torus are defined by the 
vanishing of the specific enthalpy: ho(x±) = 0. In the limit of a narrow torus, this reduces to 
h e = (2q — 3)f2^(ft 2 —x 2 )/2 (—6 < x < b) (Goldreich et al. 1986). The Lagrangean condition 
h(x±) = also in the perturbed state reduces to the two-point boundary condition (in the 
equatorial plane) (Goldreich et al. 1986) 

= D t h = -iah + uh e x . (6) 

The second equation in (3) gives h = ia<p — $ + ik~ 1 2B(fi x , so that (Goldreich et al. 1986) 

k(a 2 (f) + ia$) + (2Ba + kh%)<j> x = 0. (7) 



In the absence of a potential $, the stability of the torus is described in terms of a 
critical rotation index for each azimuthal quantum number m. The boundary conditions 
h e = become 

ea 2 + 2Ba + kh e x = 0, (8) 

where e = k<p/ <p x . The linearized equation obtained in the limit of small a in (8) corresponds 
to the slender torus approximation (b << a) of Papaloizou-Pringle, or shallow water wave 
limit (kb « 1). About uo = 0, this obtains a critical rotation index q = y/3 for all m 
(Papaloizou & Pringle 1984). Fig. 2 shows a complete stability diagram of critical values 
q c {b/a) for for azimuthal quantum numbers m = 1,2,3.. and < b/a < 1, obtained by 
eliminating A in (8) and solving for q associated with double zeros of u by continuation 
methods (e.g., Keller (1987)). 

The presented perturbations are buckling modes, associated with the same sign of the 
radial velocity at the inner and the outer face. In contrast, two-dimensional incompressible, 
vortical modes are defined by aAip = k{2B) x ip in terms of the streamfunction ip {u = ip y 
and v = ip x ). They are generally singular with divergent azimuthal velocities when ip ^ at 
the turning point a = 0, though of finite net azimuthal momentum (if) remains continuous). 
These modes fall outside the scope of the present discussion. 

A quadrupole buckling mode emits gravitational radiation at twice the frequency uj = 
2(VL a + oo') ~ 2VL a (see Fig. 3). It develops out of an internal flow of energy and angular 
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Fig. 2. — Diagram showing the neutral stability curves for buckling modes in a torus of 
incompressible fluid, as an extension of the Papaloizou-Pringle instability to large ratios of 
minor-to- major radius b/a. Curves of critical rotation index q c are labeled with azimuthal 
quantum numbers m — 1, 2, .., where instability sets in above and stability sets in below. Of 
particular interest is the range q < 2, where the m = mode is Rayleigh stable. For q — 2, 
the torus is unstable for b/a less than 0.7385 (m — 1), 0.3225 (m = 2) and, asymptotically, 
for b/a ~ 0.56/m (m > 3) 
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Fig. 3. — The frequency diagram of the pair of waves in a buckling mode on the neutral 
stability curves of a torus of incompressible fluid. The waves on the outer face are prograde 
(positive a curves; labeled for each azimuthal quantum number m = 1,2, •••,5), whereas 
the waves on the inner face are retrograde (negative a). The dot-dashed lines refer to the 
frequency uj' as seen in the corotating frame with angular velocity Vt a = M 1 / 2 /a 3 / 2 of the 
torus at major radius r = a, where the highest (lowest) curve refers to m — 1 (m — 5). 
Note that up to b/a ~ 0.3, uj' remains close to zero. Hence, the observed frequency of the 
gravitational radiation as seen at infinity is close to mfl a for low m. 
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momentum from the inner to the outer face of the torus, in which total energy and angular 
momentum of the wave remain zero. The emitted gravitational radiation is therefore not ex- 
tracted from the kinetic energy in this pair of waves. This contrasts with radiation extracted 
from single waves of frequency < uo < mQ T by the Chandrasekhar-Friedman-Schutz in- 
stability. It may be noted that this is equivalent to a positive entropy condition SS > 
in the first law of thermodynamics — SE = £It(—5J) + T5S for a torus at temperature T, 
upon radiation of waves with specific angular momentum 5J/5E = m/u to infinity. See also 
Schutz (1980) on the entropy condition in the Sommerfeld radiation condition. We proceed 
with the effect of gravitational radiation-reaction forces as follows. 



4. Gravitational radiation-reaction force 

The backreaction of gravitational radiation consists of dynamical self-interactions and 
radiation-reaction forces (Thorne 1969; Chandrasekhar & Esposity 1970; Schutz 1980). For 
slow motion sources with weak internal gravity, the latter can be modeled by the Burke- 
Thorne potential as it arises in the l\ post-Newtonian approximation 

<$>bt = \xjx k (V - l -I8^ , (9) 

where 1^ = f T.XjXkdxdy denotes the second moment tensor of the matter surface density S. 
This intermediate order does not introduce a change in the continuity equation (such as in the 
2-nd order post-Newtonian approximation (Chandrasekhar & Esposity 1970; Schutz 1980; 
Schutz & Verdauger 1983)). In cylindrical coordinates (r, 0) (x = rcos0 and y = rsin0), 
and for harmonic perturbations rj = r\(^ l6 ~ lu)t = rje 2ie e~ lU}t of the wave-amplitude, we have, 
in the approximation of a constant surface density S, 

I XiXj = S / / XiXjdxdy. (10) 

JO Ja+x-+r)- 

Explicitly, I xx = (7rE/2) [(a + x + ) 3 r] + — (a + x_) 3 ri_] , which determines = {x 2 I xx + 
y 2 I yy + 2xyI xy )/5 = z 2 I xx /5, where z = x + iy. Here, I xx = I xx e~ lbJt and z 2 comprises e 2%e 
- combined, $ = ^ e 2te -' luJt . The harmonic time-dependence e~ luJt derives from the integral 
boundaries in (10), and hence applies to all components of the moment of inertia tensor. The 
linearized radiation-reaction force derives from the fifth time-derivative, i.e., $ = — iuo^^bt 
in the stability analysis of the previous section, supplemented with the kinematic surface 
conditions — iai] = <p x on the inner and outer boundaries x — x±. Explicitly, we have 
icr ± $(x±) = i(3(l + x±) 2 K(x-, x+), where K(x-, x+) = (1 + x + ) 3 (p x (x + ) - (1 + x^) 3 (f) x (x^) 
and (3 = 7ra£(wa) 5 /10. A fiducial value (3 ~ 10 -4 obtains for for a torus mass m T = 0.1M 
and a radius a = 3M. Figure 4 shows the destabilizing effect of (3 = 2 x 10~ 4 . 
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Fig. 4. — The complex frequency diagram of the frequency oo' of the quadrupole moment in 
the torus in response to the radiation-reaction force. The results are shown for a canonical 
value of b/a — 0.2 and (3 = 2 x 1CT 4 , corresponding to a torus mass of about one percent 
of the black hole mass. The dot-dashed curves are the asymptotes for (3 = 0. This shows 
that gravitational radiation backreaction-forces renders contribute to an unstable quadrupole 
Papaloizou-Pringle buckling mode. 
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The luminosity in quadrupole gravitational radiation, C gw ~ (32/b)(VL T M) ll l z (5m T /M) 2 , 
results from a mass inhomogeneity bm^jM of a few permille (van Putten 2001). The 
Papaloizou-Pringle instability may form a quadrupole shock on the inner and other face 
of the torus. Shocks amplitudes rj± will saturate by dissipation oc r)±. Balance with energy 
input in the suspended accretion state, being on the order of the luminosity in gravitational 
radiation oc r/±, leaves rj± self-consistently to be on the order of dmr/mr of a few permille 
as well for canonical values m T = 0.1 M and M = 1OM . 



The index q is related to the pressure P in the torus. For a three-dimensional density 
p, the symmetric case in forementioned slender torus limit h = Pj p = (2q — 3)Q 2 (b 2 — x 2 )/2 
gives 



where kT =5MeV denotes a typical temperature of the torus in suspended accretion. Hence, 
q deviates from 3/2 by order unity. Furthermore, t\t — fnl{ a + fw + /h)> where fn,fw 
denote fractions of magnetic flux in the inner and outer torus magnetospheres and a = 
2qz — 2f 2 (van Putten 2001). Here, z represents viscosity between the inner and outer face 
of the torus in terms of the square of the ratio of radial turbulent magnetic flux to net 
poloidal magnetic flux. Fig. 2 shows z = 0.56, derived from the number of unstable wave 
modes in the limit of a slender torus. With a symmetric flux distribution /# = f w = 1/2, 
we have Ht/^h — l/4g. At a critical Rayleigh value q — 2, the torus emits quadrupole 
radiation at 600Hz for a 7M Q black hole. 



Black hole-torus systems may form a class of transient sources with emissions in several 
channels, possibly in association with GRBs. For long GRBs, this suggests an event rate of 
about 1/yr within a distance of lOOMpc in association with supernovae. 

Continuous all-sky supernovae surveys may detect all or most nearby hypernovae. Their 
coordinates provide a guide to searches for the proposed bursts of gravitational radiation, 
by LIGO/VIRGO, bar (e.g. (Poggi 2002)) or sphere detectors. This approach reduces the 
geometrical search space, and may be augmented by surveys for transient sources in other 
wave-lenghts, e.g., in the radio (Levinson et al. 2002). 



5. 



Estimate of index of rotation 
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6. LIGO/VIRGO searches for bursts in hypernovae 
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Searches for gravitational waves may find enhanced positive detection probability as- 
sociated with supernovae, as the sum of false positive P(+\— ) and true positive detections 
P(+|+), and in excess of P(+| — ) away from supernova events. If all long GRBs stem from 
hypernovae, we expect a local supernovae-to-hypernovae ratio of a few hundred. The excess 
P(+|+) becomes apparent for a signal-to-noise ratio better than two. Alternatively, a priori 
selected hypernova candidates from the parent supernova sample, e.g., by detection of late 
time radio afterglows, would give an apparent excess P(+|+) at lower signal-to- noise ratios. 
Complemented with distance estimates to supernovae or hypernovae, this provides a test 
for Kerr black holes (van Putten 2001). We may test for GRBs from rotating black holes 
by comparing the durations of bursts of gravitational radiation with with the de-redshifted 
durations of long GRBs. 
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